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Abstract. - In this paper we present a novel approach combining linear response theory 
(Kubo) for the conductance and the Density Matrix Renormalization Group (DMRG). The 
system considered is one-dimensional and consists of non-interacting tight binding leads coupled 
to an interacting nanostructure via weak links. Electrons are treated as spinless fermions and 
two different correlation functions are used to evaluate the conductance. 

Exact diagonalization calculations in the non-interacting limit serve as a benchmark for our 
combined Kubo and DMRG approach in this limit. Including both weak and strong interaction 
we present DMRG results for an extended nanostructure consisting of seven sites. For the 
strongly interacting structure a simple explanation of the position of the resonances is given in 
terms of hard-core particles moving freely on a lattice of reduced size. 



Introduction. During the past decade improved experimental techniques have made 
production of and measurements on one-dimensional systems possible [1], and hence led to an 
increasing theoretical interest in these systems. Since its formulation in 1992 [2] the Density 
Matrix Renormalization Group method (DMRG) has been established as a very powerful, 
quasi-exact method for numerical calculations of properties of (quasi) one-dimensional sys- 
tems. 

In this paper we present a new approach for calculating linear response conductance for one- 
dimensional interacting nanostructures coupled to non-interacting tight binding leads. The 
method combines Kubo expressions for the conductance with numerical DMRG calculations 
and is valid for arbitrary interaction strength. It facilitates a unified description of strong and 
weak interactions and provides conductance directly form a transport calculation, without 
relying on relations between equilibrium and transport properties. 

We employ current-density and current-current correlation functions to calculate the con- 
ductance and in the non- interacting case compare to exact diagonalization calculations. 

In the strongly interacting limit a simple interpretation of the position of the resonances 
is given in terms of freely moving hard-core particles on a reduced size lattice [10], and 
quantitative comparison with numerical DMRG results shows good agreement. 
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Fig. 1 - One-dimensional interacting nanostructure with Ms sites, coupled to non-interacting tight 
binding leads. The total system size is denoted M, the number of lead sites is Ml ■ The interdot and 
interlead hopping elements are trj t and t respectively, while the contact between the nanostructure 
and leads are via it and tn. 



Model. We are interested in studying the effect of correlations on transport within 
a microscopic model of an interacting one-dimensional nanostructure coupled to two non- 
interacting tight binding leads, as shown in Fig.^ Electrons are treated as spinless and only 
nearest neighbor interaction is considered. The corresponding Hamiltonian is 

Ho = i?NS + Hi, + He, (1) 

TL2 — 1 712 — 1 

hns = Y,u g c] Cj + (- WsVi + 4-is) + y sV]-is-i)' ( 2 ) 

j=ni j=ni+l 

m-1 M 

= -{^(4^-1+4-1^)-* E ( C k-l+4-lCi), (3) 

i— 2 2=712 + 1 

= — ^i( c ni £ V»i— 1 C 7ii-l C 7ii) ~~ tn( C n2 C n 2 -l + C n 2 -l C n 2 ) 

+TVV'(4 1 c ni cJ ll _ 1 c ni _ 1 + c n 2 c n 2 c i 2 -i c n2-i)- ( 4 ) 

The parameter 7y controls the smoothing of the interaction on the dot over the contact links 
as discussed in [3], and U g is a gate voltage on the structure. In this work we set t = ^Dot = 1- 

Kubo Expressions. - Using linear response in applied source-drain voltage, VsD{t), the 
current is given by 

H = H + SH, (5) 
{J n {t)) = J-if dt'(^ Q \[J n (t),SH(t')}\ij ), (6) 

J — DC 

Jn(t) = -itn[cl(t)S n -l(t)-ci-l(t)Zn(t)], (7) 

where Ho is the Hamiltonian in eq. ^ the applied voltage perturbation is 5H(t) = Vso(t)N, 
A(t) — e tHot Ae~ tHot denotes the interaction picture time evolution of the operator A, and 
\ipo) denotes the ground state. Note that in this approach A{t) contains all correlations of 
the unbiased structure, since we apply it to the quasi-exact ground state given by the DMRG 
procedure. The number operator is taken as a symmetric combination of the left and right 
lead operators, N = \(N L — Nr), and J is the equilibrium current included for completeness 
and henceforth neglected in all numerical calculations. 

2 . ~. 

The Kubo conductance in the DC limit, g = %-(«/) /Vsd, can be expressed in terms of two 
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different correlators, 

n> (H - E Q y + if 

e 2 , r Stt^Hq - E Q ) f , . 

5./ J = -r(^o Km— — -aJn 2 Vo)> (9) 

n [(i? -S ) 2 + ?7 2 ] 

where the positions rt\ and 77,2 arc defined in Fig. ^ Analogous Kubo expressions were used 
by Louis and Gros in [4], where Quantum Monte Carlo calculations for the density-density 
correlator were performed. 

Damped Boundary Conditions. - To improve the finite size scaling and to facilitate the 
use of sufficiently short leads we use exponentially damped boundary conditions, decreasing 
the hopping elements towards the end of the leads exponentially as shown for the right lead 
in cq.[ini( 1 ) 

\-t, -t, -*] -> [~t, -td, -td 2 , • • • , -tcL** ' 1 , -td M °] , (fO) 

where d < 1 . The improvement of the finite size scaling relies on two properties of the DBC's: 

(1) They allow for use of a smaller r\ and (2) serve as a particle bath for the nanostructure. 
The first property is caused by the introduction of exponentially small energyscales in the 
system thus reducing the finite size level splitting at the Fermi energy at half filling. The 
second property can be understood from the fact that the energy cost of adding or removing a 
particle from the damped region is of the order of the exponentially small hopping element. ( 2 ) 

The DBC's introduce two more parameters in the model, the number of damped bonds 
Mo and the damping factor d, and these must take values such that physical quantities do 
not depend sensitively on the particular choice. 

Numerical Calculations. Before actual numerical calculations can be performed the 
parameters of the model, Md, d, and r/, must be determined. This is done using exact 
diagonalization calculations for the non-interacting systems, specifically the resonant value 
at U g = 0.( 3 ) For fixed Md we do indeed find a range of d values that produce essentially 
identical physical results, indicating the range of validity of the DBC's. Additionally we find 
that the actual value of Md is not significant (for reasonably large values) as long as the 
corresponding value of d is tuned such that the damping at the edge reaches values of the 
same order of magnitude. The leads used are sufficiently long to keep the damped region 
separated from the nanostructure, thus allowing Friedcl oscillations at the structure edge to 
decay before reaching the damped region. 

The magnitude of the parameter rj is bounded by physical arguments; from below by the 
fact that it should be larger than the finite size level splitting to allow transport, and from 

( x )Modified BC's in connection with DMRG were introduced by Vekic and White in [5] using soft boundary 
conditions to reduce finite size effects. Note that exponential damping corresponds to the hopping Hamiltonian 
in the Numerical Renormalization Group, which models the logarithmic discretization. 

( 2 ) In principle properties (1) and (2) of the DBC's could be obtained by using longer non-damped leads. 
However these leads would have to be exponentially long making such a direct approach impossible. 

( 3 ) Considering structures consisting of an odd number of sites has the advantage that the central resonance 
(by symmetry) remains at U g = for half filled leads. Due to the bath property of the DBC's it is safe to 
assume that half filling is maintained in the parts of the leads that are close to the nanostructure. In contrast 
the strongly damped regions act like particle baths and therefore cannot maintain half filling for non-zero 
external potential. 
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above by the broadening of physical results by any finite 77, and should thus be much smaller 
than the width of the resonances we wish to resolve, ft is important to note, that r\ is an 
inherent property of any transport calculation and can only be avoided if one finds a way to 
obtain transport properties from equilibrium properties. 

The conductances in cqs. |H1 and are given in terms of ground state correlators and 
hence DMRG is directly applicable. To evaluate the correlators we use the correction vector 
DMRG [6-8] in the zero frequency limit. Calculating, e.g., the correlator in eq. OH is done by 
formulating the linear problems, 

7- -Jn^o) = |&) Jn^o) = [Ho-Eo+iriW^j), (11) 

H - E Q + If] 

which can be solved for \4>j) by a linear solver. Having solved for the correction vector \<fij) 
the conductance is found as the vector overlap, 

|^) = \<t>f) + i\4>]), (12) 

gjj = ~{4i\^)- (is) 

In our DMRG calculations we target apart from the ground state also the real and imaginary 
parts of the two correction vectors, \4>i) and \4>2)i as well as the states N\i/)o) and J ni 2 |'0o) 
to ensure that the DMRG basis is suitable for describing the conductance accurately [7,8]. 

It should be mentioned that the damped boundary conditions make the convergence rate in 
numerical calculations much slower. In addition any finite external gate voltage, U g , changes 
the particle number in the structure and the excess particles come from the bath property 
of the DBC's. We therefore face the problem that the damping should be sufficiently strong 
to provide a reasonable particle bath but at the same time a strong damping decreases the 
coupling of the highly damped region to the rest of the system. To remedy the slow convergence 
in the DMRG calculations we turn on the damping in steps and perform several finite system 
DMRG sweeps for each such damping step. In other words, we perform a complete finite 
lattice calculation employing typically 11 sweeps and then initiate the scaling sweeps. This 
allows DMRG to gradually optimize the basis to include the damping in the leads and provides 
a more gradual decoupling of the damped regions from the rest of the system, thus improving 
the convergence rate at the cost of more DMRG iterations. 

Nevertheless the resolvent equations, eq. El arc still ill-conditioned and standard solvers 
like the Conjugate Gradient Method do not converge. We use instead a preconditioned David- 
son type solver similar to Ramasesha [9] modified with a Gauss-Seidel enhanced block diagonal 
prcconditioner. The DMRG calculations presented in Fig. [2] were done using up to m = 1200 
states. In our DMRG implementation we do not fix the number of states per block to be 
m but rather fix the dimension of the target space to be at least m? . In the calculations 
presented this corresponds to an increase of block states of typically 15% — 30%. 

Results. - Here we present DMRG and (in the non-interacting limit) exact diagonal- 
ization calculations for a single resonant level, Fig. |2(a)| and a nanostructure consisting of 
seven sites coupled symmetrically to two non-interacting leads. For the extended structure 
we present results in the non interacting limit, Fig. 2(b)| and for weak and strong interaction, 
Fig. [2(c]1 and[2(dy 

The spinless single resonant level is generically non-interacting and serves as a testing 
ground for the approach. The exact result for the conductance in the symmetrically coupled 
case can be shown to be a Lorentzian of full width 4i' 2 at half maximum, where t' = = t^. 
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(a) Single resonant level, Mg = 1 and M = 
102. f's denote exact diagonalization results, g's 
denote DMRG results, and L denotes the exact 
Lorentzian in the infinite lead limit. The inset 
shows an enlargement of the resonance peak. 




(b) Seven site nanostructure, Ms = 7 and M = 
150, in the non interacting limit, V = 0.0. Exact 
diagonalization calculation. 
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(c) Seven site nanostructure, Afg = 7 and Af 
150, in the Luttinger Liquid regime, V = 1.0. 



(d) Seven site nanostructure, Mg = 7 and M = 
150, in the charge density wave regime, V = 5.0. 



Fig. 2 - Conductance, g, and number of particles on the dot, Nd, versus external potential U g for 
a single resonant level and for an extended nanostructure consisting of seven sites, gjj denotes the 
current-current correlator, and gjjN denotes the current-density correlator. The left/right contact 
hopping elements are t L / R = 0.5 and the parameters of the DBC's are Mb = 30 and d — 0.8. For 
the calculations above we use r\ = 1/M. For the interacting spectra notice the offset of resonance 
positions of the order V as compared to the noninteracting case. The interaction on the nanostructure 
is smoothed over the contacts with yv = 0.5. 



In Fig. 2(a) we show exact diagonalization and DMRG calculations for the single resonant 
level and the two sets are virtually indistinguishable. This verifies that the truncation error 
introduced by the DMRG is negligible. Furthermore we have plotted the exact Lorentzian 
result, and the agreement between the three curves is very good, demonstrating the accuracy 
of our combined Kubo and DMRG approach. 

There is a systematic difference between the current-current and the current-density cor- 
relators, specifically close to resonances the current-density correlator generally gives better 
results. This is due to the additional energy dependent broadening given by Hq — Eq in 
the current-current correlator. The opposite is true in the tails where the current-current 
correlator is more reliable since it is less sensitive to changes of the particle number. 

In numerical calculations the parameter 77 is always finite making the expected form of 
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the conductance peaks that of an area normalized Lorcntzian (La) of half-width r\ convoluted 
with the "bare" physical result. Assuming as a first approximation that the latter is a height 
normalized Lorcntzian (Lh) of width T, the expression for the expected numerical results is 
of the general form 

(14) 

To leading order in rj the conductance at the resonant level is then given by grcs 2ry/r, 
which demonstrates that one needs small r\ to reach the unitary limit, gres = 1- However 77 is 
known from the input and T can be extracted from the results. Thus a conductance value on 
resonance of 1 — 2?y/r is explained entirely by the broadening by the finite leads and therefore 
suggests that infinite leads in this case would yield the unitary limit. Our calculations indicate, 
that the peak width is only slightly affected by the interaction on the nanostructure, as long 
as the nanostructure remains in the Luttinger liquid regime. However, once the structure is 
driven into the charge density wave regime the peak width decreases rapidly. A more detailed 
study of the resonance shapes is considered future work. 

The position of the resonances can be described by the addition spectrum, 

UNd-I^Nj, = E No-l^ E N D ^ (15) 

where Eq d is the energy of the isolated nanostructure occupied by Nd particles. In the large 
interaction limit the kinetic energy of the particles can be approximated by freely moving 
fermions on an effective lattice of size Mg = Ms — Nd- In this approximation one describes 
the interacting fermions by effective hard-core particles of the size of the interaction range, 
compare [10,11]. Thus cq. EDcan be expressed as 

N D No-1 

n— 1 n— 1 v 7 

where Ad should be small enough that the nanostructure is still in a delocalizcd state. 

In an effective charging model the additional splitting of the levels due to the interaction 
is linear in the charging interaction V. By contrast, in our microscopic model the interaction 
leads to an overall offset for the non-central peaks, while their mutual splitting is governed 
by the kinetic energy, ~ t. 

In tab. H| we show a comparison of resonance positions as predicted by the reduced lattice 
(RL) model in cq.^J as predicted by exact diagonalization (ED) of the isolated nanostructure, 
and resonances found in our DMRG calculations for interaction strengths V = 5, 20, 30. The 
position of the outermost resonance from — ► 1 particle fits fairly well for both predictions, 
while the next ones deviate somewhat. The RL prediction for the transition 2 — > 3 is not 
expected to be accurate since Nd = 3 is a localized charge density wave like state. All exact 
diagonalization predictions are correct to lowest order in t/V as expected. 

Conclusion. - In this work we have presented a new approach for linear conductance 
calculations of interacting one-dimensional nanostructures, combining linear response for con- 
ductance and DMRG. We have benchmarked this new approach against exact diagonalization 
calculations in the non-interacting case and found excellent agreement, which serves as a real 
test for the real space DMRG. For the resonant level we also compared our results to the exact 
Lorentzian result, and found excellent agreement. 

For the interacting case we have presented conductance curves for a seven site nanostruc- 
ture in both the Luttinger Liquid (V = 1) and the charge density wave (V = 5) regimes, 



(L A *L H )(x) = t;t~ 
2 (x 



-7 + I72 

x y + ( v + r/2) 2 
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Table I - Table of peak positions for the Ms — 7 site structure with interaction V = 5, 20, 30, as pre- 
dicted by the reduced lattice (RL) model, by exact diagonalization (ED) of the isolated nano structure, 
and as found from the conductance peaks in our DMRG calculations. The RL prediction for Nd = 3 
is not expected to be accurate since the nanostructure is in a localized charge density wave like state. 
Except for the RL prediction for V = 5, No ~ 3, all predictions are correct to linear order in t/V . 
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thus demonstrating the versatility of our approach. We find the largest conductance when the 
particle number in the structure fluctuates, in agreement with physical intuition. 

In the large interaction limit we have shown that a simple picture based on effective hard- 
core particles moving freely on a reduced size lattice describes the position of the resonances 
quite well. However, the peak width is strongly decreased by strong interaction. 

We expect that further finetuning of the method and numerical parameters will lead to 
significantly more precise results facilitating calculations for more complicated structures and 
allow to quantitatively describe resonance peaks for strongly interacting and extended struc- 
tures. 
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